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PREFACE 



This paper presents the results of a. study of imposing a three- 
dimensional perturbation on a barotropic westerly current. The 
purpose of the present investigation was to present theoretical 
solutions to the following bs.sic problems: (1) the determination 

of the three dimensional structure of the wave; (2) the determination 
of the speed of propagation of the wave; (3) the value of considering 
vertical motion in the wave motion. 

Undertaken as the thesis requirement for the degree of Master of 
Science in Aerology, this paper was prepared at the United States Haval 
Postgraduate School, Monterey, California during academic year 1951 - 1952 . 

The author is particularly indebted to Associate Professor F. L. 
Martin of the Department of Aerology, for his advice and guidance during 
the entire preparation of this study. 
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I. INTRODUCTION 



The theoretical study of the "behavior of large-scale flow patterns 
in the extratropical zones of the earth, associated with migratory long 
waves in the "belt of prevailing westerlies, occupies an important place 
in meteorological literature. 

Many informative theoretical and synoptic studies of the behavior 
of these large-scale flow patterns have "been made since V. Bjerknes [5]. 
in 1916, advanced a theory, based on general hydrodynamic considerations, 
that cyclones originate as dynamically unstable wavelike disturbances in 
the westerly current. In 1923. H. H. Clayton made the first 

systematic use of zonal index by means of the trend of mean pressure 
versus latitude to determine increasing and decreasing zonal circulation. 
By 1933. V. Bjerknes and collaborators £27} had computed the zonal 
velocity distribution from zonal temperature and pressure distribution 
at the ground, showing the maximum horizontal velocity normally occurs 
just below the tropopause and concluded that there were 2-3 waves around 
the 60° latitude circle and slightly above 4 around the JO 0 circle. In 
1937. J. Bjerknes |j3,28] offered a simple explanation for the displace- 
ment of perturbations superimposed on zonal distribution. Variations 
in latitude and pressure perturbations, with simultaneous transport of 
air, result in changes of amplitude or longitudinal displacement of 
long waves. This result was obtained by combining the equations of 
motion with the equation of continuity. Rossby m pointed out two 
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basic objections, first, the omission of acceleration terms which 



may be of the same order of magnitude as the centrifugal and de- 
flecting force, and second, the displacement of observed flow patterns 
caused by isallobaric systems resulting from divergence associated with 
initial wind distribution. Rossby C 27 ] extended, the J. Bjerknes paper 
reasoning that since the depth of the atmosphere was much smaller than 
the horizontal scale of large disturbances, it was permissible to treat 
the atmospheric motion as two-dimensional, non-divergent flow. Further, 
he deduced that between the longer westerly moving waves, there existed 
an intermediate wave length, such that its perturbation remained 
stationary. The speed of propagation of moving waves is dependent on 
the strength of the westerlies and the wave length, and was based on 
the consideration of a change in vorticity in a vertical air column 
which has been displaced from one latitude to another. Haurwitz [l6, 17] 
extended Rossby 's 1939 solution by studying the effects of lateral limits 
of the wave and a variable Coriolis parameter on a curved earth. 

Ilamias and Clapp C223 investigated Rossby' s wave speed formula 
statistically. Their verification was generally successful and they 
concluded that amplitude and zonal velocity were both important in 
determining the stationary wave length. 

V. Bjerknes had, by this time, found the number of waves around the 
h5° latitude circle to be between four and six. In 19^» J. 3jerknes 
and Holmboe Cl9;i derived a formula. a,nalagous to that of Rossby for the 
more general barotropic atmosphere which agreed qualitatively with that 
of J. Bjerknes' of 193? • 
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By 19^7, Gharney [5] in his baroclinic atmosphere theory, based, 
on J. Bjerknes and Holmboe 1 s theory of wave motion in a baroclinic 
atmosphere, introduced a solution to the barotropic wave, ,<ith this 
theory, he obtained a solution, based on the principle that a wave 
will travel with such a speed that the pressure tendencies arise from 
the displacement of the pressure pattern in accordance with the field 
of horizontal divergence. Charney found this physical explanation 
acceptable, but the dynamics of the wave could not be solved by 
analytic methods based on semi-empirical considerations of gradient 
wind. Gharney integrated, the equations of motion, assuming 

motion to he adiabatic, and thus determined the speed of propagation, 
stability criteria and the three-dimensional wave structure. In 19^8, 
31iassen [12] generalized the Rossby-Holmboe wave formula, hut he 
neglected vertical acceleration and the Coriolis parameter when it 
appeared in multiplicative combination v/ith the vertical velocity. 

In 19^9 » Charney and Sliassen C7] extended numerical analysis by use 
of the quasi-geo strophic approximation. 

In a field where horizontal divergence was so small that the 
absolute vorticity of air parcels v/as conserved and the scale of 
motions was large enough that the Goriolis parameter variation with 
latitude v/as important in determining wave dimension, Cressman C 9 I 
scrutinized the Rossby wave velocity formula and found retrogression 
in a series of quasi-stationary long waves began when the basic zonal 
velocity decreased in speed or shifted south. Such retrogression in- 
creased the long wave number (i.e., formed a new cold trough). 
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Continuing his previous investigation, Oharney [?. 8i found, that 
the mean flow with respect to pressure could he described as two- 
dimensional non-divergent in an equivalent barotropic atmosphere* 

He used this method as a basis for partially successful twenty-four 
hour forecasts of the 5^0 mb flow pattern by numerical integration 
of the non-linear vorticity equation* The simplified equations 
derived expressed conservation of the vertical component of absolute 
potential vorticity. The hydrostatic equation together with the 
individual rate of change of vorticity evaluated geostrophically t 
reduced to a single partial differential equation for pressure 
tendency* The basis of correspondence between the barotropic and 
baroclinic atmosphere lay in the notion of an equivalent barotropic 
atmosphere, in which horizontal motion at a. particular level approx- 
imated the actual motion* 

It can readily be seen that a great amount of effort has been 
placed on the long wave theory in an effort to gain a forecasting 
tool* The use of vertical motions has suffered in con^arison, but 
primarily because of the sparsity of observations and the necessity 
to rely upon indirect information* Investigations of the vertical 
component of motion were undertaken as far back as 1911* when V* 
Bjerknes [13] described a method for such computations. In 1913* 
Bxner suggested the high correlation between 9 km pressure and 
temperature could be accounted for by vertical motions. Hesselberg, 



w 



in 1915 , extended Fxner's argument “by showing that the presence of 
persistent vertical motions of 1-10 centimeters per second would 
account for the observed magnitude of pressure-temperature changes. 
Jeffreys in his 1922 study of vertical velocities necessary to 

mainta-in hailstones, found vertical velocities of any magnitude at 
all occurred only in eddies and could be ignored in treating mean 
motion. In 1927, Haurwitz showed that in an atmosphere of 

constant lapse rate, that the maximum vertical component of velocity 
was in the region of maximum divergence and was of the order of eight 
centimeters per second. Brunt, [4} in his Memoirs Q-93^1 found vertical 
velocities involving divergence were of the order of a f ew centimeters 
per second. 

As observational information increased, quantitative analysis of 
divergence and vertical currents became more practicable. jjurst M. 
in 19^0, found vertical velocities of a few centimeters per second; 
Fleagle [xj] t ia 19^5 » Gaviola and Fuertes C15J, in 19^7 » assumed 
adiabatic processes and computed the vertical component of velocity. 
Panofsky [237 • 19^6, from third degree polynomials, found values 

of the same magnitude (i.e., a few centimeters per second). 

.ihile vertical motions of several meters per second have been 
observed, their duration was short and smoothed out over a twelve - 
hour period. Thus, the high valued measurements by theodolites as 
observed by Suring M have not altered the basic contention that 
persistent vertical motions are of a much smaller order. 
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Oharney T6], in assigning magnitudes to vertical motions in 19^8, 
explained that the failure of the tendency equation in pressure prog- 
nostication was due to the discounting of vertical notion. Gharney £3], 
further detei’mined that forecasts with the quasi -geos trophic assumption, 
could he achieved only when the distribution of vertical motion and sea- 
level pressure distribution we re known. Thus, it is necessary to extend 
numerical forecasts to permit forecasting vertical motion. In 19^9» 
Bellamy £2^ described a method of objective calculation of vertical 
velocities. In 1951> Phillips [25] indicated that the equivalent 
barotropic atmosphere was incapable of providing distribution of vertical 
motion and sea-level pressure adequately for forecasting the "weather n . 

In view of the interest in vertical velocity distribution, it is the 
purpose of the present investigation to present a theoretical solution 
of a two-layer barotropic, incompressible atmosphere with three-dimensional 
velocity divergence zero, introducing vertical motion as a perturbation. 
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II. FORMULA T 1 017 ..ID RESULTS OF TH3 PROBLEM 

A brief summary of the principal subject matter of the investigation 
unhampered by mathematical detail is now given to set forth the procedure 
followed and results obtained. 

Chapter III concerns the construction of a barotropic model. The 
troposphere is characterized by a constant vertica.l lapse rate; the 
stratosphere is assumed to be isotherma.1, with constant zonal wind 
throughout the atmosphere. The equations of motion are given in ordinary 
cartesian coordinates for an incompressible atmosphere. Horizontal 
velocity divergence is assumed to be compensated by vertical motion. 
Boundary conditions are formulated. It is shown that the model is 
consistent with the boundary conditions and the steady state notions 
prescribed. 

The last section of Chapter III, treating the actual motion as a 
small perturbation superimposed on an undisturbed constant zonal current, 
is now discussed. The velocity components, as well as the pressure 
component, are expressed a.s simple harmonic perturbations of infinite 
lateral extent traveling in a west-east direction at constant velocity. 
These perturbation components are functions of height as well as of 
horizontal coordinates and time. The vertical velocity, involving a 
general expression of the z-coordina.te, is used as a basis for the 
determination of all other velocity components. To simplify the 
mathematical difficulties, the meridional perturbation amplitude is 
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made proportional to that' of the zonal perturbation. The zonal 
amplitude, in turn, hy the prerequisite of zero of velocity divergence, 
is simply related to the vertical motion amplitude function, H(z). It 
is shown that pressure is also a function of the vertical motion ampli- 
tude function. 

Since the vertical amplitude is in terms of H(z), the amplitudes of 
the other components are determinable from H(z). The determination of 
the proportionality factor between the meridional and latitudinal per- 
turbation velocities assumes an important role in the solution of the 
problem. Upon its determination rests the dependency of the wave 
velocity and the maximum allowable meridional perturbation velocity 
for a particular maximum zonal perturbation velocity. By analogy to 
Rossby’s wave velocity, a restriction that the proportionality factor 
be less than unity is assumed. It is shown that the proportionality 
factor is dependent upon the discriminant of an auxiliary quadratic 
equation. Following the boundary conditions this discriminant must 
have a value slightly in excess of unity. The discriminant value 
must be determined from an inequality. The choice of the value of 
this inequality plays a dominant role in any further solution of 
this problem. 

The problem then becomes one of determining the. speed of propa- 
gation of the wave and the character of the vertical perturbation 
which introduces the amplitude as a function of height. The deter- 
mination of these are dependent upon wave length and the parameters 
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characterizing the mean state of the atmosphere, namely, the height 
of the tropopause, the vertical velocity at the tropopause, the 
zonal wind, the vertical lapse rate, the mean latitude and surface 
density, pressure, and temperatiire. 

In Chapter IV, the form of the perturbation is set forth. The 
perturbation functions are solved in the stratosphere, and the tropo- 
sphere (by approximation) . The proportionality constant between 
meridional and zonal perturbation velocities is determined. This 
constant is of great importance in the determination of the wave 
speed and stationary wave length. Amplitude evaluation is determined. 
The complete solution of the approximate case is shown diagramatically. 

In Chapter V, the tropospheric perturbation functions are solved 
without simplifying the coefficients of the differential equation 
satisfied by the amplitude function, H(z). 

Due to barotropy, no phase change is introduced into the structure 
of the wave. Vertical cells do appear in the perturbation velocity 
structure moving with the basic wave velocity. The perturbation 
fields of velocity diminish with increasing height and eventually 
approach zero. The wave in the meridional velocity field lags 90° 
behind the wave in the pressure field and the wave in the density field 
is in phase with the pressure wave. The zonal perturbation is in phase 
with the pressure wave, whereas the vertical velocity perturbation leads 
the pressure wave by 90°. These phase relationships are shown in 
Figure 3. 
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III. TEL. MuTHSKiiT ICAL F ORlfJL^TIOL 

1. The Atmospheric Model. 

As an approximation to the atmosphere in middle latitudes, a 
model of the undisturbed state is characterized as follows: (a) the 

undisturbed wind is zonal; the speed is constant and independent of 
the coordinates throughout the troposphere and stratosphere; (b) the 
lapse rate of temperature is constant in the troposphere and zero in 
the stratosphere; (c) both troposphere and stratosphere are barotropic, 
incompressible and frictionless layers. 

2. Fundamental Equations. 

It is assumed for purposes of mathematical simplicity that the 
curvature of the earth can be neglected. This assumption is va.lid when 
length of wave is small compared with circumference of the zonal circle 
along which the wave moves. This planar motion is expressed in a 
rectangular system of coordinates x, y and z with x increasing to the 
east, y northward, and z vertically upward. Then u, v, w denote 
corresponding total velocity components (consisting of undisturbed 
velocity and perturbation velocity). Further, denote total pressure 
by p, angular velocity of the earth by.fi, acceleration of gravity 
(assumed constant) as g, geographical latitude by ^ , the z and y 

components of Coriolis acceleration, f andy^ respectively. The 
density , having been assumed incompressible, is constant and 
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equal to the unperturbed value . he jJulnrLan equations of motion 
in ordinary cartesian coordinates, fixed in the earth, become 



sLl. s - J- Sfc t f ~/% '> r , 
dt ^ > 



dor - -_L - f 

dt /« 



( 1 ) 



( 2 ) 



duS - -d~ ~d)P + /3 CL 

d-t ' ^ 9Z ' 
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A fourth eqjuation is obtained from the assumption that the three- 
dimensional velocity divergence is zero. 



^ *V r o . 






A relationship between density and height in the troposphere (with 
the assumption that g is constant with height) provides the fifth 
equation 
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But the density-height relationship in the stratosphere is 
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where h^ is height of tropo pause. 
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3. The Boundary Conditions 



The bo\u>dary conditions express the following physics.! properties 
of the motion: (a) the perturbation components are continuous across 

the tropopause, that is, 

Ap r A u. = Av - A Cur = A p - AT - O. 

(b) the momentum vanishes at the limit of the atmosphere. 

|im /0<Z - Itm / oV' =■ 1 1 km CZr = O . 

z->*> / ■*.-*•» ' •*-**» ' 

k. Derivation of the Perturbation liquation. 

The fundamental equations (1-6) together with the boundary con- 
ditions (7), (8) impose the necessary restrictions on the theoretical 
model, .vith the convention of a capital letter for steady sta.te value 
and a small letter for the same quantity in the perturbed state, the 
characterization of the steady sta.te is given by the conditions 

W = dU = ^ EE - o , U = constant. (9) 

“3T s>t 9x 

Equation (2) implies that the undisturbed flow must satisfy the 
condition of geostrophic equilibrium 

o = -i-"3 P-.fl/. (10) 

"ay 

Equation (10) implies that yo and P are functions of y and z only. 
Equation (3) expresses the condition for hydrostatic equilibrium in 
the mean state 
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with, the omission of \J , which can be shown both empiric? lly and 
theoretically to be negligible in comparison to the acceleration of 
gravity. 

By means of the equation of state, and equations (4) and (5) with 
prescribed values of U, P 0 T 0 , T, and h f|1 , the theoretical steady state 
model is completely described. The boundary conaitions(7) and (8) are 
a.lso satisfied by our steady state model. Hence the model prescribed 
is consistent with the fundamental equations (1-6) and the boundary 
conditions (7), (8). 

5. The Linearized equations. i 

The motion to be investigated can be treated as a small perturbation 
with velocity components u, v, and w superimposed on an undisturbed 
constant zonal current U, with the assumption of infinite lateral 
extent. Thus the motion is described by 

u = \i(x,z,t) + U, v = v(x,z,t), vJ s w(x,z,t) . (12) 

Similar expressions for the density and pressure in the disturbed state 
are 

? = p(x,z,t)+ P(y», jo. ^(y, z) . (13) 

In deriving the perturbation equations the usual assumptions are 
ma.de. (1) The tota.l disturbed plus undisturbed motion satisfies the 
hydrodynamic equations, as well as the undisturbed motion alone. 

(2) The perturbations are so small that terms of second order in the 
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perturbation quantities can be neglected with respect to terms of 
the first order in the perturbation quantities. Linearizing the 
'iulerian equations of motion (1-3) and simplifying by means of the 
steady state equations (9-11), the system (1-3) becomes: 

^ u, 4. U + A kjS - f'O's-l ~d Jg , (1*0 

9* 9x 7 

f O jhc + -P « O 

3£ 3* 

B UJ~ -f (J — /£ U- : -i_ g) )g. 

^ 3 X / ^ 9 Z 

including the velocity divergence, 

9 * 



4 jlLdT - O . 
9Z 
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IV. THE BAROTROPIC PERTURBATION 



1. The Form of the Perturbation and Resulting Conditions. 

The perturbation components, traveling in the x- direction at a 
constant velocity, are assumed to be simple harmonic motions of infinite 
lateral extent. The amplitudes of the perturbation are included in a 
general expression of the z-coordinate, as follows 



u = sin m(x-ct) F(z); v ■> cos m(x-ct) Q(z); 
w = -m cos m(x-ct) H(z); p = sin m(x-ct) C(z). 



(18) 



To simplify the mathematical difficulties, the assumption is 
made that F(z) is proportional to Q(z), D being the constant of 
proportionality. Since u and w of eouations (18) are related by the 
condition of zero velocity divergence, it can be shown that 

F(z) = H 1 (z ) = Q(z) D. (19) 

The equation for the velocity of the wave is derived by substituting 
the perturbation velocities (18) into the equation ( 15 ) • iith the 
equalities of ( 19 ) known, the wave velocity equation is 

c = U - fD (20) 

m 

Since U and m are real, c will be either real or complex, depending 
on D. The determination of D will be discussed fully in Chapter IV-2. 
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Substituting the perturbed, velocity, pressure and mass fields into 
the perturbation equations (16) and (14) leads to 



^ -y> -f O m H C-z) + Hci) (21) 

^ PfYOjifj Ufa +/yS^(z) y (22) 

respectively. 

The differentiation of (22) with respect to z, leads to 





By differentiation of density from equation (5) , we have 



i ■ 

Equating (21) and (23) and replacing X> and by their 

92 

expressions in (5) and (2^) leads to the following ordinary 
differential equation for the amplitude function H(z): 

( I- o u)H(z) -K Hbt) - j2u£Ci-<<5]«^>= 9 (25) 

where K is an abridged notation for 

K C2/fir-i) ■ 

Using the further abridged notations, 

A- AKVrn. 8 = 

■f C v^-\ > v 



(16) 



equation (25) > the tropospheric differential equation becomes 

(l-u z) H^z ) - K H'fr) iSA t H(z> =<=> * ( 28 ) 

2. Perturbation Functions in the Stratosphere, 

In exactly the same manner as described in the preceding section, 
the stratospheric differential equation follows by using (6) instead 
of (5) for the density expression. Thus, analogously to equation (2^), 
the stratospheric density derivative is ' 




where h^ = RT ^/g is the height of a homogeneous layer at strs-tospheric 
temperature T^. s/e, therefore, obtain the second order linear differ- 
ential equation with constant coefficients 

The auxiliary equation of (30) is 



hZ 



dW: 

D-) 



(3D 



whence 



% = 



_ j_ 

2-hi 



r - 4 -ts, I hfyPmh*- ~( 2 Pjzitjtl . ( 32 ) 

- *VV TU^ i r < d 0 



The determination of the value of D is dependent upon the value 
assigned to the discriminant of the auxiliary equation (32). By 
analogy with the Ho s shy wave velocity 



U 






(i?) 



it would follow that, for Ross hy waves, 



0 R - /3/Bfm ^ o.iv • 



The only restriction placed on D s,t this stage of the investigation * 
is that 1, from which it follows at once, that the discriminant of 
(28) has a value n^ , at latitude ^5°* slightly in excess of unity, 

nj>l. 

Solving the discriminant for D: 



-6 £- l')"]-£* l 02J) X»t -l (33) 



Solving equation (33) by the quadratic formula gives 



Y* +y* (1 




-'T' 




. Yx J X ll 1 


j • 



In order to expand the radical in (35) "by means of the "binomial 
expansion, the condition that 



must hold. 

* Dyl v/ith a moderately long wave of 6000 km, would give a value 
of U-c > 100 mps, so that obviously D ? 1 does not apply. 



let 




(34) 
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Examining the two possible cases: (a) 1, (t>) Y^7 1. 

Case (a): Y^^l. 

3 quation ( 36 ) becomes 

(2 m hj’' ? (n i -\) /nl j 

so that 

I ?1 < \ , ( 36 ') 

3or L = 6000 km, < 1.012, Then from the inequality of 
case (a), the expansion of the radical of ( 35 )i v/ ith the choice of 
the negative sign, yields 



D = (Vfr-l) *- -t-J- . o?) 

8 YJ 

Case (b): Y^> 1 
Equation (3^) becomes 

{2fnhX/(nl-l) y 1 . 

It follows that 

J < K?z < t 1 ^6 M ^ ^ (38) 

where n^< 1.0004 for a wave length of 6000 km. The value of 
determined by case (h) is included in that of case (a), so that Y ^<.1 
will he used in the evaluation of D. Moreover, all terms of (37) beyond 
the second will he considered negligible. Thus at latitude 45° where 

^ ^ 
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( 39 ) 






The evaluation of D , which for a particular mean latitude and wave 
number remains constant, is dependent on the value of n 2 selected 
from the inequality ( 36 ' )• In this investigation, it was assumed 
desirable to have the ratio of the absolute maximum velocity of the 
horizontal perturbation to that of the vertical perturbation as large 
as consistently possible. This ratio, hereafter called 



M - I > (39’) 

I urmoK \ m 

can be shown to be a function of n ^ and L. The value of D is com- 
pared in Table 1. for various values of and wave number 1J. 



m in 
meters - -*- 



H 


n a. 


x 10 4, 


2 rn h ^ 


*2. 


M 


D 


4 


1.0056 


0.888 


0.0113 


0.107 


0.5 


.376 


5 


1.0070 


1.110 


0.0141 


0.119 


0.5 


• 

“O 

00 


6 


1.0084 


1.332 


0.0170 


0.132 


o.5 


.376 


7 


1.0098 


1.55^ 


0.0198 


0.141 


0.5 


.377 


8 


1.0112 


1.776 


0.0226 


0.151 


o.5 


.377 


9 


1.0125 


1.998 


0.025^ 


0.161 


o.5 


.377 



Constant parameters: = 45°, T z = 218° A 

Table 1. 
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Having determined the value of D qualitatively in terms of m, 



and n , the general solution of the stratospheric equation (30) is 

+C^ Z (t0) 

with the abridged notation 

- '/zhz. • 

By application of the "boundary condition (8), yOuT— it follows 

Z-*<* J 

that M<z) — » o if and. only if the coefficient C is made zero, 

2->»© 2 
The particular solution satisfying the "boundary conditions is 

y(z) =: C / e S ^ l ~" z) ^ . (41) 

In the determination of the pressure perturbation in the stratosphere, 
the amplitude function G(z) is rewritten from equations (22) and (6), in 
the abridged notation, as 



Gd) ft fO-gJO-nJ -y,) 

From the solution of (4l), the perturbation components of the 
stratosphere are defined. In abridged notation, they are 

U z - C i S^Chn^) Slh 

- C»_ S^O-O e s ^ /_/7 ^ z Cos m(x-ct) (43) 

p 

ur z ^ m C, e ^ 1 h ^ 7 ~ cos mCx-c-t) ? 

- Q (z) si n Yn(x-c$> m 



(21) 



In order to solve (43) for Cq, suitable numerical values must be 
assigned to the constant parameters , hj, T g , w^ y , H and n 4 , /° . 
The following values are selected: ^ = 45°, h,j = 10 km, T^ 218°A, 
w m ax = 5*5 cm sec~l, H = 5. n x = 1.00?, /° 0 - 0.001121 gm cm“3. 

With assignment of these quantities, numerical values may be derived 
for the parameters m, Sj_, M, where m = 1.11 x 10"^ meters - !, 

Sa. = .0?85 km - l, M - -g-. liquation (43) yields Cq = -5 x 10^ m*" sec - ^. 

3. An Approxima-te Solution in the Troposphere. 

To obtain the solution of the barotropic wave by approximating 
the value of (1 -*z) in the tropospheric equation, (28) let 

(j ■ xZ ) — /'To (44) 

be a satisfactory approximation. Substituting this approximation 
into equation (25), the equation becomes a second order differential 
equation with constant coefficients, namely. 



% H 'Yz\ - K H'Cz) -h 

~rf 

r =(!-£) $//?> 0 < £ < | 






UCz) = o. 



(45) 

(46) 



Then, if h 3 = RT 3 



equation (45) can then “be shown to he 



H'(z) - € H'cz) f- / /3 £ jjjzl— 
^3 h 3 £p 2 --i] 



- D*n^ Wr}~o. (47) 
D 1 - 1 J 



( 22 ) 



The general solution of this equation is dependent upon the discriminant 
obtained from the auxiliary equation of equation (47). By letting the 



discriminant equal , where 




( 48 ) 



with D known from the stratospheric solution^ Hj is determined 
Then the general solution in the approximate case is 




(49) 



4. Application of the Tropopause Boundary Condition, 

The method of evaluating the coefficients C^, and is the same 
for both the approximate and actual case. The approximate case is no w 
considered. 

'.(ith the boundary condition of a discontinuity of the first order 
at the tropopause the perturbation components must be continuous across 
the tropopause, 




*5 • 



(50) 



Then from (49) the first derivative with respect to z is 




( 23 ) 



so that 




(52) 



From the boundary condition (50) , equate vertical velocities and also 
horizontal velocities at the tropopause, thus 



By solving these two equations simultaneously C 2 . ana C^ are determined. 
5. The Complete Solution in the Approximate Case. 

A r ith C^ and C^ determined, the perturbations of the troposphere 
are by approximation defined as 

^3 ' u 3(ma\) Siw nn ( J 




(53) 

C, S t C a S t Ci-n^)i r= S 3 e 5 ’ >>r[C 3 c / + qG -%) « V] 



o 



u 3<m«x; Cos ™ Ct-ct); 



*^3 " Scmox) ^ os 



(5^) 




(24) 



In order to solve (5^)* suitable numerical values must be 
assigned to the parameters of the troposphere, in addition to those 
assigned in the stratosphere (e.g. , Section 2 of this Chapter). The 
additional parameters, having selected values, are 
Tj = 256. 5 0 A,f* = 5.18° km' 1 , dith assignment of these quantities it 
can be shown that the constants C a and have the values 3*59 and 
-4.98 x lcA, respectively. A diagram of the perturbation velocities 
is presented in Figure 1. for a graphical representation of the 
solution determined in this Chapter, 
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(26) 



V. TEE VARIABLE COEFFICIENT CASS 



1. The Solution of the Perturbation. 

with D determined and constant throughout the atmosphere for 
any particular mean latitude, the second order differential equation (28) 
is now solvable. By suitable manipulation, described in the Appendix, 
this equation (28) is a case of the confluent hyper geometric differ- 
ential equation 



y L"(p) Cb~?l L'(p) - L(p) - o , 



(55) 




Equation (55) is satisfied by the functions 




f - (yj‘' b F(q.~\3 + i j Z'b, y) - 



I + ~ b fi 




The confluent hypergeometric equation then has the general 



solution 




( 57 ) 
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To test the convergence of ^ — F*Ca.^b y ip)' » compare to the 
series Q * , where \ y\ • it can he shown that, termwise, 

f P f < and orders higher than the second can he dropped with 

accuracy to the first five significant figures. Then 

% * i +-%-»? +- ~p h 2 - (53) 

' h r bCbbi) ' j 

where a - 

Testing the convergence of 

(-p)^ b ^ y 

after the first 12 terms, the series is more rapidly convergent than 
where X* - % Ipl. In the range of values used in this study, however, 
there is negligible error, if series is shortened to 




I + 0.- bfj ^ j_ (ci-bll) fa ~b Wp t? 3 - 

2 (A-6)(3-b) 



(59) 



The general solution in the troposphere from equation (57) is 

HU) . + 

C q c-p) , ~ b f l-h &Tb+t y +. j.(x-b-H%l~b + z) h 7] . ( 60 ) 

x 2 (aCTJWbT- ' J 
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2. Application of the Tropopause Boundary Conditions. 

After the manner of evaluating coefficients in Chapter IV, 
equating the vertical and horizontal velocity across the tropopause, 
two resultant equations can he solved simultaneously for C z and C^ . 

3. The Complete Solution. 

,/ith Cj_and C^ determined, the perturbations in the troposphere 
are defined as 



L4 X --(wi tc } [?i(-r)'' h +(-!>ryo- *») ^j) s» ma-cti 

=. -w (C L f J f- C 3 Cos 

' V X ~ + i fz y'l 1 ’' 1 *] ^i)j Cc>s ">(X-C i)j 

D m x 



(61) 



(-?r b /«-b) % ]) +/3 (c t f, +c 3 vjl s'» ^c^ L i 



yp ~ D m ^ Dm _ j~/ _ or 



( 29 ) 



In (6l) the primes denote differentiation with respect to z. The 
values of parameters to he inserted in (6l) are 
D = .377, m - 1.1 x 10“'-' cm - -*-, ot = 1.9 x 10" cm“^, 

K = 1.064 x 10" cm - - 1 -, so that p can he evaluated for any value of z. 
In particular, at the tropopause p has the value p = -.01?6. Like- 
wise p - .415 x 10 - °, and is constant throughout the troposphere. 

The complete solution includes the solution above the tropopause 
which remains precisely as developed in Chapter IY-2. 
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VI. Discussion OF RESULTS 



It would appear that the variations of the value of np determined 
"by the discriminant of the stratospheric auxiliary equation ( 32 ) plays 
the most important role in this entire investigation. By variation of 
this parameter alone, the wave velocity, the stationary wave length 
and the ratio of horizontal to the vertical maximum velocity, can he 
varied, within the allowable range of this parameter, wave velocities 
determined by Rossby are included. 

Gravity and sound waves are also included in these solutions for 
appropriate values of n. The values of n correspond to the inequality 
: l tTrh x ./L- > with values of the wave-length, L, 
suitably chosen to be representative values of wave-length for sound 
and gravity waves, respectively. It will be remembered that Charney 
£5] , in dealing with the problem of the unstable baroclinic wave, 
first eliminated these short stable waves, which in the present paper 
is equivalent to restricting the condition to D < 1. 

A further possible investigation is suggested by this study. By 
limiting the lateral extent of the disturbance, larger values of the 
velocity components are consistent with the small values of vertical 
divergence used in this study, the values of and 3 a 7 ' 

would then partially compensate one another, giving the proper (small) 
magnitude to < ^h' ^ - 2bif . 
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aPPSITDIX I. 



The transformation of equation (28) to the form of a confluent 
hypergeometric differential equation follows, 
aquation (28) is restated here. 



)H(Z) -K HCz.) + (A + eO-*zj J WU) =o. 



( 28 ) 



Divide "by (A+B): 




« 

MB 




- o. 



Divide '*(t-ecv/ A + B )^0 J 



then multiply "by (A+- 3) 



(28. a.) 



I - OtZ 

1 -JBgcz. 
A+B 



H(Z) 




«&> * (Aw) = o. 



(28. b) 



Momentarily, consider the first term only 







(28. c) 



(35) 
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♦ I • . 
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Expanding, 




if it is assumed that higher order terms negligible. 
Equation (28. h) together with (28.c,d.) yield 




Collecting terms, (28. e) is 



Hfa- +- M H<i> »°- 

Let 



5 - i g<2. 

A + A+5 

then 



l~g dfj = -A* 
A ' dZ A+B 



J -/- _£*z 

At3 

and equation (28. f) is 



I 4- 



A 



BJ- 

4 



* 



- K(| +-!■!) -♦ (A4 e) wet) - o , 



(28. d) 



(28.e) 



(28. f) 



(28. g) 



( 36 ) 



where 



Hoo = *]j = J(D J 

hU) - AL±L ~ - jH)A* 

d<$ Az A+ 8 

iquation (28. g), in terms of £ , becomes 

(me)^ 



- O. 



(28. h) 



Dividing (28. h) by +&J yields 



By letting r = KBCA+BjAu 
equation (28. i) becomes 

§J"(f)+ v (fi+£_ - S) JCjp)+0fB/7c§)=o. (2B,J) 

N B y 73 30- 



( 3 ?) 
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Then if \T £ - fp > 

^ ?UL) L (>?) } 

T<J) r &sf - " 

n?> = r -V(#-)' *"L"(p)' 

equation (28. j) is 

1 - 5)l'c ? ) i Mb£ l O) --o. 

0 (« /l) 1 ' 

Divide equation (28. k) by /* and substitute the equality | 

the equation becomes 

^/V +fi"CA±Bl -/! L ' f ^ LM = o. 

a J v - f */*) 1 ' 

if b* r*(A+&)/e> , a = -fa+B?/ r *( e (Aj L » 

equation (28.1) is identically equation (55) » which is the form of 
the confluent hypergeometric differential equation. 



(28. k) 
, and 

( 28 . 1 ) 



( 38 ) 
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